Spacetime noncommutativity and ultra-high energy cosmic ray experiments 
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If new physics were capable to push the neutrino-nucleon inelastic cross section three orders of 
magnitude beyond the standard-model (SM) prediction, then ultra-high energy (UHE) neutrinos 
would have already been observed at neutrino observatories. We use such a constraint to reveal 
information on the scale of noncommutativity (NC) Anc in noncommutative gauge field theories 
(NCGFT) where neutrinos possess a tree-level coupling to photons in a generation-independent 
manner. In the energy range of interest (10 10 to 10 11 GeV) the ^-expansion (\6\ ~ I/A N c) an d 
therefore the perturbative expansion in terms of Anc retains no longer its meaningful character, 
forcing us to resort to those NC field-theoretical frameworks involving the full #-resummation. 
Our numerical analysis of the contribution to the process coming from the photon exchange, pins 
impeccably down a lower bound on Anc to be as high as around up to 900 (450) TeV, depending 
on the estimates for the cosmogenic neutrino flux. If, on the other hand, one considers a surprising 
recent result occurred in Pierre Auger Observatory (PAO) data, that UHE cosmic rays are mainly 
composed of highly- ionized Fe nuclei, then our bounds get weaker, due to the diminished cosmic 
neutrino flux. Nevertheless, we show that even for the very high fraction of heavy nuclei in primary 
UHE cosmic rays, our method may still yield remarkable bounds on Anc, typically always above 
200 TeV. Albeit, in this case one encounters a maximal value for the Fe fraction from which any 
useful information on Anc can be drawn, delimiting thus the applicability of our method. 

PACS numbers: ll.10.Nx; 13.15.+g; 13.60.Hb; 98.70.Sa 



The observation of ultra-high energy (UHE) neutrinos 
from extraterrestrial sources would open a new window 
to look to the cosmos, as such neutrinos may easily escape 
very dense material backgrounds around local astrophys- 
ical objects, giving thereby information on regions that 
are otherwise hidden to any other means of exploration. 
In addition, neutrinos are not deflected on their way to 
the earth by various magnetic fields, pointing thus back 
to the direction of distant UHE cosmic-ray source can- 
didates. This could also help resolving the underlying 
acceleration in astrophysical sources. 

Since the GZK-structure in the energy spectrum of 
UHE cosmic rays at ~ 4 x 10 19 cV has been observed 
recently with high statistical accuracy [l|, the flux of 
the so-called cosmogenic neutrinos, arising from photo- 
pion production on the cosmic microwave background 
PIcmb — > A* — > Nir and subsequent pion decay, is now 
guaranteed to exist. Although estimates for the cosmo- 
genic neutrino flux are very model-dependent, primar- 
ily due to our insufficient knowledge of the nature and 
the origin of UHE cosmic rays, possible ranges for the 
size of the flux of cosmogenic neutrinos can be obtained 
from separate analysis of the data from various large- 
scale observatories [2], b| ■ With this knowledge of cosmo- 
genic neutrinos and the increased exposure of large-scale 
detectors, a tight modcl-indcpcndcnt constraint on the 
neutrino-nucleon cross section can be placed |4[ . A simi- 
lar analysis but with a different approach to the problem 
was carried out in [5[. 

We note that the uncertainty in the flux of cosmogenic 
neutrinos may be further aggravated in view of appar- 
ently contradictory recent results regarding the chemical 
composition of UHE cosmic rays. While the Pierre Auger 



Observatory (PAO) data indicate that UHE cosmic rays 
may be highly- ionized Fe nuclei [6J, the results from the 
High Resolution Fly's Eye Collaboration (HiRes) are con- 
sistent with the proton dominance in the UHE region 
[7[ . In the former case the flux of cosmogenic neutrinos 
would possibly be much smaller than that in the latter 
case [8(. Further confusion was brought in by the fact 
that both experiments are intrinsically inconsistent re- 
garding the arrival directions of UHE cosmic rays [9[. 
Note also that the most recent evidence of anisotropy of 
UHE cosmic-ray arrival directions at the 99% confidence 
level [10| by the PAO would question further its inter- 
pretation that primary UHE cosmic rays are likely to be 
dominated by heavy nuclei. On top of that, there is an- 
other inconsistency besetting the PAO data: the average 
of the depth in the atmosphere at which the shower con- 
tains the largest number of high energy particles implies 
an Fe fraction < 60%, whereas the analysis of its stan- 
dard deviation implies for the same fraction to be > 90%. 
Although previously most observations and arguments 
seemed to favor protons as a dominant component at 
the highest energies, and although some speculations ex- 
ist that inconsistencies in the PAO data rather point to 
some new physics than to a heavy component in the UHE 



cosmic rays 
121. 



Ill ] , this issue still remains an open question 



In this paper, we employ the aforementioned upper 
bound ,4j on the cross section to constrain the scale of 
noncommutativity entering NCGFT in which neutrinos 
couple directly to photons. We use upper bounds de- 
rived by non-observation of UHE neutrino- induced events 
in RICE Collaboration search results [13[ and exploit- 
ing the predictions for the cosmogenic neutrino flux from 



[2j, |3j. Also, to realize the second possibility as indicated 
by the PAO data, we derive our bounds by using the 
same predictions for the cosmogenic neutrino flux, but 
now diminished by the fraction of FE nuclei in the UHE 
portion of cosmic ray spectrum, which we treat here as a 
free parameter. 

NCGFT that we consider is based on Moyal-Wcyl 
star(*)-product [4l| and it belongs to a class of mod- 
els that expand the NC action in 9 before quantiza- 
tion [lj,ll5(- In these models Scibcrg-Wittcn (SW) maps 
[ 1 €31 ] are necessary to express the noncommutative fields 
ip and A^, that appear in the action and transform under 
NC gauge transformations, in terms of their asymptotic 
commutative counterparts ip and A p . The algebra gener- 
ated here is the enveloping algebra. Computation of any 
(physical) quantities, obtained in this approach, are given 
as a power series in the NC parameter 9, thus showing 
no UV/IR mixing on the quantum level. 

Admittedly, in other approaches, with #-unexpanded, 
i.e., models based on ^-products only [17[, the celebrated 
quantum-gravity effect shows up: the UV/IR mixing [18l — 
|20| . The same would happen in our approach. 

Coupling of matter fields to Abelian gauge bosons is a 
NC analogue of the usual minimal coupling scheme. It 
will imply only a vector-like NC coupling between pho- 
tons and neutrinos in the U(l)* gauge-invariant way, thus 
endowing also otherwise sterile right-handed (RH) neu- 
trino components a new interaction. 

It is possible to extend the model [15J to the NC elec- 
troweak model based on the other gauge groups, like NC 
version of the standard model (NCSM), etc. 



21, m 



That NCSM [21( at (9-order does break Lorcntz symme- 
try, appears to be anomaly free [23| | and it does have 
remarkable well behaved one-loop quantum corrections 
1241 , as well as some other models do I25I427II . 



Signatures of noncommutativity and/or the bounds on 
the NC scale come from neutrino astrophysics la, [28J, |2£ 
and cosmology [3C 



31 



from high energy particle physics 
[32H34 1 , producing a scale of noncommutativity of order 
few TeV's. Typical low energy nonaccclcrator experi- 
ments are the Lamb-shift [35[ , and clock-comparison ex- 
periments [36| . The bound from Z — > 77 decay on Anc 1 
of order a few TeV's 331, is the most robust due to the 



finite one-loop quantum corrections in the gauge sector 
of the nmNCSM [24j. Also, the interesting upper bound 
on the scale of NC, of order of one tenth of the Planck 
mass, can be obtained from the black-hole physics in the 
NC spacctime 37 1. 

Neutrinos do not carry an electromagnetic charge and 
hence do not couple directly to photons. However, in 
the presence of spacetime NC, it is possible to couple 
neutral particles to photons via a *-covariant derivative 
D^ip. Thus, the action for neutrinos that couples to an 
Abelian gauge boson in the NC background, invariant 



under the NC gauge transformations, is 

S N c = d?x($*ij*D M $-m$*i>) , (1) 

D$ = d^-iKelAf,*^}, (2) 

Ap = An + eerPApidvAp-ldM + OiP), (3) 



ip = v + 






°A p d u ij) + 0{9 2 



(4) 



where Q and Q are, up to the first order in 9, ordi- 
nary SW maps for the Abelian NC gauge potential and 
"crural" SW map for the NC "crural" neutrino field, re- 
spectively [42]. Here, we fix the NC scale by defining 
qpv = c^^/A^fQ, such that the matrix elements of c are 
of order one. 

Performing the SW and ^-product expansions up to the 
first order in 9 one obtains the following electromagneti- 
cally gauge-invariant NC action for photons and neutri- 
nos in terms of commutative fields and parameters [la ] : 



S 



3 = -- J <t l x$A llv (i9>""'d p -0>"'m)il>, (5) 



with 9^ vp = e» v Y + Q vp l^ + Q p ^l u \ ^V = d^Av-dvAp. 
From ([5]) we extract the Feynman rule for massless 
left/right(L/R)-handcd neutrinos (the same for each gen- 
eration and k = k' + q) [l5| : 

T h {Dvi) = l e ( 1± ^ r ^' (6) 

The double differential deep inelastic scattering cross 
section d 2 o~nc I dxdy for the process vN — > v + anything 
via the vv^f vertex as given by (j6]) , can be calculated with 
the help of the structure functions, 
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((coi - cia) 2 + (co 2 - c 23 ) 2 ) -fr^ xy (1 - y), (8) 
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dxdy E u Mj\r(xy) 2 

< [(1 - V)F2 + y 2 xF? + y(l - y/2)xF2] 



(9) 



Here the Bjorkcn variable x is given in terms of the inelas- 
ticity y, the total energy of the system s and the invariant 
momentum transfer —Q 2 as x = Q 2 /(ys). The azimuth 
angle ip in ([5} is defined with respect to the direction of 
the incident neutrino. The quantity <7 comm refers to the 
cross section a neutrino would have as if it were a Weyl 
particle with a charge e. The structure function F^s for 
the isoscalar nucleon are given by 



2xF?(x,Q 2 ) - F^{x,Q 2 ) 



(10) 



^(x,Q 2 ) = \(Ql + Q 2 d )[u{x,Q 2 ) + u(x,Q 2 ) 

+ d(x,Q 2 )+d(x,Q 2 )} 

+ 2Q 2 c c{x 1 Q 2 )+2Q 2 s s{x 1 Q 2 ) 

+ 2Q 2 b b(x 7 Q 2 )+2Q 2 t(x,Q 2 ) 



Fj(x,Q 2 ) = 



(11) 
(12) 



and in our analysis we use the parton distribution func- 
tions from the audibly available package CTEQ4 j38| . 
The kinematical region of interest here is that of high Q 2 
and very small x values x w 1.7 x 1Q~ 7 /{Ev/10 11 GcV), 
because of the rapid increase of the parton densities to- 
wards the small x's. Different extrapolation approaches 
can result in uncertainties as large as a factor of two at 
E V = IQ 12 GeV |H. 

Employing the upper bound on the vN cross section 



derived from the RICE Collaboration search results [13 1 
at E v = 10 11 GcV (4 x 1(T 3 mb for the FKRT neu- 
trino flux [2| ) ) , one can infer from (|oT ll2[) [also employing 
(coi — C13) 2 + (C02 — C23) 2 = 1] on the scale of noncommu- 
tativity Anc to be greater than 455 TeV, a really strong 
bound. One should however be careful and suspect this 
result as it has been obtained from the conjecture that 
the 0-expansion stays well-defined in the kinematical re- 
gion of interest, and the more reliable limits on Anc are 
expected to be placed precisely by examining low-energy 
processes |2jj. Although a heuristic criterion for the va- 
lidity of the perturbative ^-expansion, v^/Anc % 1, 
with s = 2E„Mn, would underpin our result on Anc> 
a more thorough inspection on the kinematics of the 
process does reveal a more stronger energy dependence 
EJ s 1//4 /Anc ^ 1- In spite of an additional phase-space 
suppression for small x's in the ^-contribution 34[ of the 
cross section relative to the ^-contribution, we find an un- 
acceptably large ratio a(9 2 )/a(9) ~ 10 4 , at Anc = 455 
TeV. Hence, the bound on Anc obtained this way is 
incorrect, and our last resort is to modify the model 
adequately to include somehow the full-0 resummation, 
thereby allowing us to compute nonperturbativcly in 9. 
The simplest possible modeling, which we propose here 
is: (a) Replacement of the NC fields with the commuta- 
tive ones, i/> — > %j) and A^ — > A^; (b) Expansion of the 
•-product up to all orders, and resummation to the 9- 
exact final results. Following (a) and (b), the NC action 
instead of (JlJ can be written now as 



S NC (0) = -ie / d 4 ^7"(A M ^-^A„). (13) 



above action corresponds to the 

H. Next 



At this point the 

neutrino-photon interaction proposed in Ref 
we expand the ^-product to all orders in 9, and take 
derivatives on ip(x) ~ e tkx tp(k) and A^x) ~ e tqx A f _ l (q). 
After getting the momentum dependence order by or- 
der we resume the obtained series into the nontrivial 
exponential phase factors exp(±i^2~), and extract a 



generation-independent Fcynman rule for masslcss L/R 
neutrinos 



rfoM 



ie(l±7 5 )7^sin( — ) 



(14) 



By comparison of our Feynman rule (|14|) with those 
from [18|, [20( , it is clear that our model would also pro- 
duce the UV/IR mixing in computing quantum correc- 
tions, thus, our model is not perturbatively rcnormaliz- 
able and it is not clear how to intcrprcte the quantum 
corrections and to relate them to the observations [40(. 
In return, we obtain the 9 well-behaved deep inelastic 
cross section at ultra- high energies. 

With the aid of the full vertex (fl"4|) , the neutrino cross 
section still retains the form (JT]), but with the relevant 
quantity now defined as 



1 f 2lx 
X = — / d<p4sm 2 ( 

2tt Jo 



kck' 

zil NC 

2(l-cos(A)J (B)) , 



(15) 



where Jo is the Bessel function of the first kind of order 
zero, and 



A = " " co 3 (cos#- 1) 

A NC 

E V E' V . 
B = —2 — sin it 



(16) 



xsign(c i - c 3)\/(coi - C13) 2 + (c 2 - c 2 3) 2 -(17) 

The behavior of the cross section with the scale of non- 
commutativity at fixed E v = 10 10 GcV and E v = 10 11 
GcV, together with the upper bounds depending on the 
actual size of the cosmogenic neutrino flux (FKRT [2| 
and PJ [3j) as well as the total SM cross sections at these 
energies, are depicted in our Figure 1. In order to max- 
imize the effect of noncommutativity, in our numerical 
calculations we choose Coi — C13 = C02 — C23 = C03 = 1. 

A few comments are in order regarding qualitative be- 
havior of the cross section. The common attribute for 
both energies is the existence of a plateau at small Anc's, 
where the cross section tends to a constant value. Also, 
the oscillatory term in (fT5")) enters the regime of rapid os- 
cillations for lower Anc's for both energies. This behavior 
can be read off from (|15[) by noting that the argument of 
the cosine function grows as Anc is decreasing (for fixed 
energies), whilst the amplitude of oscillations as given by 
Jo tends to zero asymptotically in the same limit. The 
cross section for the upper energy does clearly surpass 
the measured value below around 900 (450) TeV for the 
FKRT (PJ) flux, bringing us a valuable information on 
the NC scale. As far as the lower energy is concerned, 
the upper bounds assigned to both fluxes are crossed, 
at energies around 350 (500) TeV, respectively. Thus, 
the conservative bound for the energy E v — 10 10 GeV is 
around 350 TeV. 



If future data confirm that UHE cosmic rays are com- 
posed mainly of Fe nuclei, as indicated, for the time be- 
ing, by the PAO data, then still valuable information on 
Anc can be obtained with our method, as seen in Fig. 
2. Here we see the intersections of our curves with the 
RICE results (cf. Fig.l) as a function of the fraction a of 
Fe nuclei in the UHE cosmic rays. We see that in all but 
one cases our method is capable to extract information 
on Anc even for Fe fractions as large as ~ 90%, giving 
always a conservative bound ~ 200 TeV. The terminal 
point on each curve gives the highest fraction of Fe be- 
yond which our method becomes unserviceable to draw 
any information on the NC scale. Thus, when more data 
on Fe fraction become available, one can easily pinpoint 
the lower bound on Anc with the aid of Fig. 2. 

To summarize, we have used the upper bound on the 
neutrino-nucleon inelastic cross section, derived by non- 
observation of UHE neutrino-induced events in extensive 
airshower arrays, to deduce information on the scale of 
noncommutativity in NC gauge field theories. The fact 
that we deal here with a process whose center-of-mass 
energy is three orders of magnitude higher than that be- 
ing achieved with terrestrial accelerators, has also im- 
plications for the theory itself. Namely, the dimension- 
less parameter involving the scale of noncommutativity 
and entering the perturbative expansion shows a strong 
incident-energy dependence, such that perturbative ex- 
pansion at such UHE energies is no longer meaningful. 
After having treated our model nonperturbativcly in 6, 
we used such a model to derive a very strong bound on 
the scale of noncommutativity, at around 900 (450) TeV, 
depending on the model for the total cosmogenic neutrino 
flux. Since the prediction for the NC neutrino-photon 
vertex in a perturbative setting is quite robust amongst 
different models, we believe that our 0-resummed inter- 
action is also a generic one. For those inclined to ac- 
cept the PAO interpretation, we have shown that even 
for the large fraction of Fe nuclei in the UHE cosmic 
rays (~ 90%) one may obtain remarkable bounds on the 
scale of noncommutativity (Anc ~ 200TcV). Finally, 
one should keep in mind that the limits derived here are 
conservative ones, meaning that any influence of Fe nu- 
clei as well as NC physics on the cosmogenic neutrino 
flux (and consequently on the experimental upper limit 
on the cross section) has not been treated. 
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FIG. 1: vN — > v + anything cross sections vs. Anc for 
E u = 10 10 GeV (red lines) and E v = 10 11 GeV (blue lines). 
FKRT and PJ lines are the upper bounds on the neutrino- 
nucleon inelastic cross section, denoting different estimates 
for the cosmogenic neutrino flux (see the text). SM denotes 
the SM total (charged current plus neutral current) neutrino- 
nucleon inelastic cross section. The vertical lines denote the 
intersections of our curves with the RICE results. 
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FIG. 2: The intersections of our curves with the RICE results 
(cf. Fig.l) as a function of the fraction of Fe nuclei in the UHE 
cosmic rays. The terminal point on each curve represents the 
highest fraction of Fe nuclei above which no useful information 
on Anc can be inferred with our method. 
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